Flooding of Chaotic Eigenstates into Regular Phase Space Islands 
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We introduce a criterion for the existence of regular states in systems with a mixed phase space. 
If this condition is not fulfilled chaotic eigenstates substantially extend into a regular island. Wave 
packets started in the chaotic sea progressively flood the island. The extent of flooding by eigenstates 
and wave packets increases logarithmically with the size of the chaotic sea and the time, respectively. 
This new effect is observed for the example of island chains with just 10 islands. 
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One of the cornerstones in the understanding of the 
structure of eigenstates in quantum systems is the semi- 
classical eigenfunction hypothesis i n the semiclassi- 
cal limit the eigenstates concentrate on those regions in 
phase space which a typical orbit explores in the long- 
time limit. For integrable systems these are the invari- 
ant tori. For ergodic dynamics the eigenstates become 
equidistributed on the energy shell 2]. Typical systems 
have a mixed phase space, where regular islands and 
chaotic regions coexist. In this case the semiclassical 
eigenfunction hypothesis implies that the eigenstates can 
be classified as being either regular or chaotic accord- 
ing to the phase-space region on which they concentrate. 
Note, that this may fail for an infinite phase space |3|. 

In this paper we study mixed systems with a compact 
phase space, but away from the semiclassical limit. Here 
the properties of eigenstates depend on the size of phase- 
space structures compared to Planck's constant ft. In 
the case of 2D maps this can be very simply stated 
a regular state with quantum number m = 0,1,... will 
concentrate on a torus enclosing an area (m + 1/2) ft, as 
can be seen in Fig.CJc). 

We will show that this WKB-type quantization rule is 
not a sufficient condition. We find a second criterion for 
the existence of a regular state on the m-th quantized 
torus, 



7m < 



(1) 



Here r H = ft/ A c h is the Heisenberg time of the chaotic sea 
with mean level spacing A c h and 7 m is the decay rate of 
the regular state m if the chaotic sea were infinite. Quan- 
tized tori violating this condition will not support regular 
states. Instead, chaotic states will flood these regions, see 
Fig.QJa). In terms of dynamics we find that wave packets 
started in the chaotic sea progressively flood the island as 
time evolves. Partial and even complete flooding is possi- 
ble, depending on system properties. These findings are 
relevant for islands surrounded by a large chaotic sea. 

We numerically demonstrate the flooding and the dis- 
appearance of regular states for the important case of 
island chains. In typical Hamiltonian systems they ap- 
pear around any regular island. On larger scales they 
are relevant for Hamiltonian ratchets p, the kicked ro- 
tor with accelerator modes Q, and the experimentally 



0, @, El and theoretically ^(J studied kicked atom sys- 
tems. The flooding of regular islands by chaotic states 
is a new quantum signature of a classically mixed phase 
space. This phenomenon shows that not only local phase- 
space structures, but also global properties of the phase 
space determine the characteristics of quantum states. 

Before we explain the origin of Eq. ([T]). we numeri- 
cally study its consequences. We choose a system, where 
we can change r H by increasing the system size without 
affecting the rates j m . A one dimensional kicked system 



H(p,x,t)=T(p)+V(x)J25(t-n), 



(2) 



has a stroboscopic time-evolution given by the mapping, 
x t +i = x t + T'(p t ), pt+i = Pt ~ V'(x t +i). The phase 
space is compact with periodic boundary conditions for 
x t G [0,M] and ^ G [—1/2,1/2]. Choosing the functions 
V'(x) and T'(jp) appropriately [Tl] we get a chain of M 
islands, one per unit cell (see Fig. [I}. The islands cover a 
relative area A 



0.215 and have fine structure close to 
their boundary that is negligible for the quantum proper- 
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FIG. 1: (a) Husimi representation of a chaotic state flooding 
the regular islands. Shown are three phase-space cells out of 
M — 1597. Full curves are classical tori close to the border 
of the regular islands and black dots are iterates of a chaotic 
orbit. This eigenstate extends well into the islands, while 
having no weight in their central region, (b) Average of the 
eigenstate over all M cells, (c) For M = 1 eigenstates concen- 
trate either on the chaotic component (left) or over the m-th 
quantized regular torus. For all plots ft e ff ~ 1/30 is used. 



ties studied here. Points inside an island are mapped one 
unit cell to the right, i.e. the island chain is transporting. 

The eigenstates of the quantum system are deter- 
mined by the eigenvalue equation, U\ip) = e^l^), where 
U is the unitary time-evolution operator over one time 
period, U = exp[— 2iriT (ft) / h e &\ exp[— 27riV(x)/h e ^]. 
The spatial periodicity after M cells requires an effec- 
tive Planck's constant h e E — M/N , with incommensu- 
rate integers M and N. We choose for M/N the rational 
approximants of /i e ff = l/(d + a) with a = (y/E — l)/2 
the golden mean and, e.g., d = 29 in Fig.[IJ This ensures 
that there are no undesired periodicities and that h e s is 
approximately constant when varying M. Moreover, the 
operator U reduces to an N x TV periodic band matrix. 
Using the symmetrized version of the map and making a 
unitary transformation to a band matrix we are able to 
calculate the eigenstates of U up to N ~ 10 5 . 

For M = 1, Fig. [IJc) snows a typical chaotic eigen- 
state and five regular states. The chaotic state extends 
over the chaotic phase-space component and the regular 
states concentrate on quantized tori. The eigenstates are 
represented on the classical phase space by the Husimi 
distribution, where for visualization we use tilted coher- 
ent states adapted to the shape of the island. For larger 
system sizes we find that chaotic states flood the islands 
of classically regular motion. Fig.JHa) shows such a state 
for M = 1597 that clearly ignores the outer tori of the 
island, which for M = 1 act as barriers for chaotic states 
(Fig. HJc), left). In the central part of the island, how- 
ever, this state has essentially no weight. This partial 
flooding of the island is observed even better in Fig.[IJb), 
where an average of the Husimi function of this state is 
taken over all M unit cells. The almost constant value 
in the chaotic component extends well into the island. 
Inside the island, clearly away from its outer boundary, 
the Husimi function sharply drops to zero. 

For a quantitative description of this flooding we now 
analyze the weight W of each eigenstate inside the is- 
lands. We determine this weight by integrating the nor- 
malized Husimi function (calculated on a 30 x 30 grid 
per unit cell) over the islands. In the semiclassical limit, 
^eff —> 0, regular states have W = 1, while chaotic states 
have W = 0. The distribution P(W) of these weights 
for all eigenstates is shown in Fig. for various system 
sizes and ft eff « 1/10 El- For M = 1 we observe, as 
expected, a main peak near W = coming from the 
chaotic eigenstates and two distinct peaks at larger W 
from the two regular states, see Fig. EJa). A remarkable 
shift of the main peak of P(W) to larger values of W 
can be observed in Fig. EJb). This shows that by en- 
larging the system size M all chaotic states continuously 
increase their weight inside the regular islands. This in- 
crease stops when the center of the main peak reaches 
W = 0.22, which corresponds to the area A reg of the is- 
land. For these system sizes all states completely flood 
the island (Fig.|2(c), inset), as observed in Ref. Q. 

What happens to the regular states as M is increased? 
Figure |5(b) shows that the corresponding peaks in the 
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FIG. 2: Distribution P(W) of the weights W of the eigen- 
states in the regular islands for /i e ff ~ 1/10. (a) M = 1: The 
main peak near W = is due to chaotic states and two fur- 
ther peaks are due to regular states m = 0, 1 (see insets), (b) 
For increasing M the main peak shifts to larger values of W 
(white line) and the two peaks from the regular states disap- 
pear sequentially, (c) M — 10946: A narrow peak remains 
around W = A re g- Three phase space cells of an eigenstate 
show a complete flooding of the islands. 



distribution P(W) disappear. Notably, the peak for 
m = is much longer visible than the peak for m = 1. 
The nemesis of the regular states can be quantified by 
determining their fraction f reg as a function of M. To 
this end we define a state to be regular when its weight 
W inside the islands exceeds 50%, where the exact cri- 
terion does not affect our analysis. Figure [3| shows that 
the fraction f reg decreases from approximately A reg all 
the way to zero. The decay is slower for smaller h e s. 

Remarkably, in Fig. [3| we see strong signatures of the 
decrease of f reg already for small system sizes M w 10. 
This holds even for small values of h e ^ where the com- 
plete flooding of the island is numerically not accessible. 
Similarly, a shift of the main peak in Fig. Efb) is clearly 
detectable for such small systems. We thus find that par- 
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FIG. 3: Fraction / reg of regular states vs. system size M for 
h e fi ~ 1/5 (triangles), 1/10 (circles), and 1/20 (squares). An 



approximately linear decrease with In M to / rf 



can be 



seen (lines are a guide to the eye). Already for small system 
sizes M « 10 and small h e f£ a significant reduction of regular 
states is observed. 



tial flooding of regular islands is easily observable. 

Why do chaotic states flood the islands of regular mo- 
tion and why do regular states disappear as the system 
size is increased? Let us consider a single regular island 
coupled by tunneling to a chaotic sea. If the chaotic 
sea is infinite, its states form a continuum. A regular 
state on the ra-th quantized torus has a decay rate j m 
to the continuum pjj. Thus, it is not an eigenstate, 
but it is dissolved into chaotic states. As a consequence, 
the chaotic states occupy the m-th quantized torus of 
the island. If the chaotic sea is finite, but large enough, 
this decay of the m-th regular state may still take place. 
The condition for the decay is that during the time l/j m 
the discrete chaotic spectrum is not resolved, leading to 
l/7m <C h/A c h = r H 0. On the other hand, if the 
chaotic sea is so small that during the time l/j m the 
chaotic spectrum is well resolved, then the regular state 
m does not decay, yielding Eq. JI}. Note, that 7 m in- 
creases monotonically with m, as for larger m the m-th 
torus is closer to the boundary of the island. 

The quantized tori of an island can thus be grouped 
into two classes: (i) The innner tori, m = 0, ...,m* — 1, 
where condition (QJ is fulfilled and regular states exist, 
(ii) The outer tori, m = m*, m max — 1, where Eq. (QJ 
is violated and which is flooded by chaotic states. Here 
m max is the number of quantized tori at a given h e &. We 
find for the fraction f veg of regular states and the weight 
W c h of chaotic states inside the island: 



freg — A 



reg 



-, W ch = A reg [l 



(3) 



Variation of the system size M in our example allows 
to change the Heisenberg time r H ~ M, while keeping 
the rates j m fixed. Enlarging M leads via Eq. to 
a decrease of m*, starting from m* = m max — 1 all the 



way to m* = 0. Together with Eq. (|3J) this explains 
Fig-Etb), where the regular state with m = 1 disappears 
before the m = state and the weight W c h grows un- 
til W c h — A Teg1 where the island is completely flooded, 
see Fig. 12c). This alsoexpl ains the decrease of f reg from 
A Yeg to 0, as observed in Fig. 01 This decrease occurs over 
an exponentially large range in M, due to the roughly ex- 
ponential dependence of j m on m. A quantitative under- 
standing requires a theory for the decay rates j m: which 
is a subject of current research on dynamical tunneling 
0,^3- Note, that in the case of chaos- assisted tunneling 
the splitting of symmetry related regular states fluctuates 
strongly, depending on individual chaotic states. In con- 
trast, the decay rate j m describes an average tunneling 
to a continuum of chaotic states. 

Variation of /i e ff affects both j m and r H in Eq. QJ. 
While r H ~ M/fteffj one expects in analogy to WKB 
theory that 7 m ~ exp[— g(m/m max )/h G s], where the 
system specific function g decreases monotonically to 
g(l) = 0. From the definition of m* follows m*/m max = 
# -1 [^eff m (^/^eff)], where g~ x decreases monotonically. 
Decreasing h e & reduces the argument of g~ x such that 
m*/m max increases. Eq. (|3J) implies that / reg grows 
and W c h decreases. Note that in the semiclassical limit, 
h e ti —> 0, we obtain m*/m max — > 1. This is in agreement 
with the semiclassical eigenfunction hypothesis, namely 
/reg = A reg and there is no flooding. In contrast, if the 
system size is infinite, we have an infinite r H and our argu- 
ment leads to m* = 0, i.e. complete flooding, for any h e fi. 
This coincides with the considerations of |3( implying a 
failure of the semiclassical eigenfunction hypothesis. 

Our explanation is complete for systems without lo- 
calization. For example, this is the case if the average 
classical drift of a unit cell is non-zero, like in atom optic 
experiments in the presence of gravity 0, ^1 • Localiza- 
tion, however, sets a lower bound to the effective mean 
level spacing, A cri ~ 1/A, where A is the localization 
length. For M > A, this leads to r H ~ A and m* stays at 
its value for M = A. According to Eq. (JHJ) the same holds 
for /reg and W c h- This applies, e.g., to dynamical local- 
ization in the kicked rotor. For transporting islands, like 
in the model studied here, A ~ I/70 is unusually large 
p3.Hl.H8]. such that already for M = A one has m* = 0, 
freg = 0, and W Q 
no consequences 

We generalize our analysis to the dynamics of wave 
packets, which is experimentally of great relevance A 
wave packet started on the m-th torus will be restricted 
to that region if condition (QJ is fulfilled, i.e. m < m*. 
If m > m*, however, the wave packet will decay into 
the chaotic sea. Particularly interesting is the case of a 
wave packet started inside the chaotic sea. The island is 
progressively flooded, i.e. the m-th torus at time t m ~ 
l/j m for m > m*. For t > r H the weight W(t) will 
saturate at W c h? Eq. (J3J). This is confirmed in Fig.^J for 
increasing values of M. 

Our results have consequences for spectral statistics in 
mixed systems which go well beyond the previously stud- 



ch = ^4-reg- In this case, localization has 
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FIG. 4: Weight W in the island vs. time t for a wave packet 
initially started as a coherent state in the chaotic sea at 
(x,p) = (0.4, -0.2) for system sizes M = 1, 13, 144, 1597 and 
h e ff ~ 1/30. The thick (green) line guides the eye to the 
linear increase with hit until W — A re g- Insets show the 
time-evolved wave packet averaged over all cells for the case 
M = 1597, demonstrating the progressive flooding of the is- 
land. 



ied effects of dynamical tunneling (see e.g. 16, 20]). The 
effective size of the regular region, f reg in Eq. J3J), enter- 
ing the Berry-Robnik formula [2l[ is drastically reduced. 

Our analysis applies as well to hierarchical states [2^ , 
which are confined by partial transport barriers with 
turnstile areas smaller than h. We predict the addi- 
tional condition 7 < l/r H for their existence, where 7 
describes the decay through these partial barriers. For 
regular states on island chains within that hierarchical 
region, condition applies, with r H given by the mean 
level spacing of the surrounding hierarchical states. 

Finally, we emphasize that the time periodicity of the 
system and the restriction of our discussion to maps 
is not crucial and that we expect flooding of islands for 
any Hamiltonian with a mixed phase space. We stress 
that this new quantum signature of chaos for eigenstates 
and wave packet dynamics already appears for small sys- 
tem sizes, e.g. island chains of length 10. This makes 
numerical explorations very feasible and should lead to 
experimental observations, for example using optical lat- 
tices 3110. 
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